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. The problem of inter-band tunneling in a semiconductor (Zener breakdown) 

in a nonstationary and homogeneous electric field is solved exactly. Using the 
exact analytical solution, the approximation based on classical trajectories is 
studied. A new mechanism of enhanced tunneling through static non-one- 
dimensional barriers is proposed in addition to well known normal tunneling 
■ solely described by a trajectory in imaginary time. Under certain conditions 

on the barrier shape and the particle energy, the probability of enhanced tun- 
| neling is not exponentially small even for non-transparent barriers, in contrast 

' to the case of normal tunneling. 
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I. INTRODUCTION 



A control of processes of quantum tunneling through potential barriers by external sig- 
nals is a part of the field called quantum control which is actively developed now, see, for 
example, Ref. [1] and references therein. Excitation of molecules, when one should excite 
only particular chemical bonds [2-4] , formation of programmable atomic wave packets [5] , a 
control of electron states in heterostructures [6,7], and a control of photocurrent in semicon- 
ductors [8], are typical examples of control by laser pulses. A control of quantum tunneling 
through potential barriers is also a matter of interest, since tunneling is a part of many 
processes in nature. The computation of probability for a classically forbidden region has 
a certain peculiarity from the mathematical stand point: there necessarily arises here the 
concept of motion in imaginary time or along a complex trajectory [9,31,34]. The famous 
semiclassical approach of Wentzel, Kramers, and Brillouin (WKB) [9] for tunneling prob- 
ability can be easily (in the case of a static potential) reformulated in terms of classical 
trajectories in complex time as a simple change of variables. The method of complex tra- 
jectories can also be applicable to a nonstationary case [12,13] which is not trivial. The 
method has been further developed in papers [14-18], when singularities of the trajectories 
in the complex plane were accounted for an arbitrary potential barrier (see also [19]). Recent 
achievements in the semiclassical theory are present in Refs. [20-24] . 

Let us focus on the main aspects of tunneling under nonstationary conditions. When the 
electric field S cos fit acts on a tunneling particle of the initial energy E, it can absorb the 
quantum Ml (with the probability proportional to the small parameter S 2 ) and tunnel after 
that in a more transparent part of the barrier with the higher energy E + Ml. The pay in 
the absorption probability may be compensated by the probability gain in tunneling. In this 
case the system tends to absorb further quanta to increase the total probability of passing 
the barrier. This mechanism of barrier penetration is called photon-assisted tunneling. If 
Ml is not big, the process of tunneling, with the simultaneous multi quanta absorption, can 
be described in a semiclassical way by the method of classical trajectories in the complex 
time [14-18]. When a tunneling particle of the energy E is acted by a short-time pulse, the 
tunneling probability is associated with the particle density carrying away in the outgoing 
wave packet. The particle energy after escape is E + SE, where the energy gain SE = Nhw, 
should be extremized with respect to the number of absorbed quanta N and the energy fko 
of each quantum [25,26]. This mechanism relates to semiclassical method for nonstationary 
potentials. 

The semiclassical method in quantum mechanics is very elegant and constructive. On 
the other hand, this method is based on the delicate mathematical issue called Stokes phe- 
nomenon when a solution is not expected to appear but it appears [27]. This makes a use 
of the semiclassical method to be not trivial even in the case of a static potential. In this 
situation the role of exactly solvable problems in quantum mechanics is very important. For 
example, in the problem of reflection of a particle from certain static potentials one can 
follow in details a formation of a semiclassical solution from the exact one [9] . Is it possible 
to find an exactly solvable problem for a nonstationary potential (excepting not interesting 
parabolic one, where the effect is trivial [17,18]) to see how the exact solution is reduced to 
complex classical trajectories under semiclassical conditions? An exactly solvable problems 
would be extremely useful since the method of classical trajectories in nonstationary prob- 



2 



lems is still challenged and an exact analytical solution would place classical trajectories in 
the rank of mathematical theorem. 

Such exactly solvable nonstationary problem exists. This is an inter-band tunneling 
(Zener breakdown) in a semiconductor [28] in a nonstationary electric field which is a con- 
stant in space. The latter condition makes the problem to be exactly solvable since in the 
momentum representation the Schrodinger equation is reduced to the first order which can 
be solved by the method of characteristics. In this paper it is studied how the exact solution 
turns over into one corresponding to classical trajectories and under which conditions this 
is possible. The goal of the paper is not to investigate Zener effect in real semiconductors 
but to use this situation to mathematically justify the method of classical trajectories. 

Another issue of this paper is that a new mechanism of enhanced tunneling through static 
non-one-dimensional barriers is proposed in addition to well known normal tunneling solely 
described by a trajectory in imaginary time. As shown in the paper, under certain conditions 
on the barrier shape and the particle energy, the probability of enhanced tunneling is not 
exponentially small even for non-transparent barriers, in contrast to the case of normal 
tunneling. 



II. PHOTON-ASSISTED TUNNELING 

A penetration of a particle through a potential barrier is forbidden in classical mechanics. 
Only due to quantum effects the probability of passing across a barrier becomes finite and 
it can be calculated on the basis of WKB approach, which is also called the semiclassical 
theory. The transition probability through the barrier, shown in Fig. 1, is 

W ~exp[-A (E)} (1) 

where 

A (E) = ^J dx^2m [V{x) - E\ (2) 

is the classical action measured in units of h~. The integration goes under the barrier between 
two classical turning points where V(x) = E. One can use the general estimate A ~ V/ftw, 
where V is the barrier height and ui is the frequency of classical oscillations in the potential 
well. A semiclassical barrier relates to a big value V/fko ^> 1 . 

What happens when the static potential barrier V(x) is acted by a weak nonstationary 
electric field £(£)? In this case there are two possibilities for barrier penetration: (i) the 
conventional tunneling, which is not affected by £ (t), shown by the dashed line in Fig. 2(a), 
and (ii) an absorption of the quantum hfl of the field £{t) and subsequent tunneling with 
the new energy E + Ml. The latter process is called photon-assisted tunneling. The total 
probability of penetration across the barrier can be schematically written as a sum of two 
probabilities 

W ~ exp -— + — ii exp (3) 



hw J \ h J \ ftw 

where £q is the Fourier component of the field £{t) and the length a is a typical barrier 
extension in space. The second term in Eq. (3) relates to photon-assisted tunneling and it is 
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a product of probabilities of two quantum mechanical processes: absorption of the quantum 
HQ and tunneling through the reduced barrier V — Ml. Since in quantum mechanics one 
should add amplitudes but not probabilities, Eq. (3) is rather schematic and serves for 
general illustration only. For example, the accurate perturbation theory starts with a linear 
£-term. When the frequency is high Q > u, the second dominates at sufficiently small 
nonstationary field (a£n/h) 2 > exp(— VL/uj). This is a feature of tunneling processes since, 
normally, a nonstationary field dominates at bigger amplitudes (aSu/h) 2 > 1. When the 
second term in Eq. (3) exceeds the first one, further orders of perturbation theory should be 
accounted which correspond to the multiple absorption, shown in Fig. 2(b). 
Let us specify a shape of a field pulse in the form 

£{t) = , (4) 

{i + t 2 /9 2 f v ; 

with the Fourier component £q ~ S9exp(— 9\Q\). In this case, in addition to the almost 
steady flux from the barrier, an outgoing wave packet (of the maximum amplitude y/W) 
is created which carries away a certain particle density. The probability of the transition 
through the process of absorption of iV quanta and subsequent tunneling with the higher 
energy E + NMl, shown in Fig. 2(b), can be estimated as 

exp [-A (E + Nhn)} = {^OaS) exp(-A) (5) 

where 

19 

A = — 5E + A (E + 5E) (6) 

Here the total energy transfer SE = NhQ is introduced (Q > 0). The maximum squared 
amplitude W of the outgoing packet is associated with a maximum value of w and is 
determined by the extreme SE to get a minimum of A. This results in the condition 
dA(E + SE)/d5E = 0. An existence of such a minimum is possible if growing up of small 
8E reduces A, that is, under the condition 29 < k \ dA (E)dE | of sufficiently short pulses. 
In other words, sufficiently short and not very small pulses (however, still much smaller than 
the static barrier field) strongly enhance tunneling by photon assistance. 

In principal, this type of treatment of nonstationary tunneling is reasonable since a 
maximum of the outgoing packet should relate to some extreme condition. Nevertheless, the 
approach (5) is schematic since a quantum interference between absorption and tunneling is 
neglected. What is a "scientific" way to account precisely a nonstationary field in tunneling? 
We go towards this in III. 



III. TRAJECTORIES IN IMAGINARY TIME 

According to Feynman [29] , when the phase of a wave function is big, it can be expressed 
through classical trajectories of the particle. But in our case there are no conventional 
trajectories since a classical motion is forbidden under a barrier. Suppose a classical particle 
to move in the region to the right of the classical turning point xt in Fig. 1 and to reach 
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this point at t — 0. Then, close to the point xt, x(t) = xt + ct 2 (c > 0) and there is no a 
barrier penetration as at all times x(t) > xt- Nevertheless, if t is formally imaginary, t = it, 



the penetration becomes possible since x(ir) 



x T 



CT 



is less then xt- Therefore, one 



can use classical trajectories in imaginary time to apply Feynman's method to tunneling. 
In the absence of a nonstationary field a classical trajectory satisfies Newton's equation in 
imaginary time 



m 



d 2 x dV(x) 



dr 2 



dx 



(7) 



where V(x) is the static barrier in Fig. 1. The classical turning point x T in Fig. 1 is reached 
at r = with the initial condition dx/dr = 0. The classical trajectory x = x(t) can be 
considered as a change of variables x — > r in the WKB exponent (2) when it becomes of the 
form 



An 



(8) 



So, in the absence of a non-stationary field, use of classical trajectories in imaginary time is 
obvious and it is simply reduced to a change of variable. 

Is it possible to extend the method of classical trajectories in tunneling to a non- 
stationary case? 

Suppose some pulse of an external field, for example (4), acts on a tunneling particle. At 
the point xt, besides an almost constant background, there is a wave packet of the outgoing 
particles shown in Fig. 3. If the method of classical trajectories is applicable, the maximum 
of the squared amplitude of the outgoing packet 



W ~ exp(-A) 

should correspond, with the exponential accuracy, to the classical action 





m 1 


if* 






~2~ \ 



dx 



+ V(x) - x£(ir) - E 



(9) 



(10) 



since a classical trajectory relates to an extreme value of the action. The classical trajectory 
should satisfy Newton's equation 



m 



d 2 x _ dV(x) 

W 2 



dx 



with the conditions 



dx 



0; 



-SUt) 



x(itq) = 



(11) 



(12) 



r=0 



We suppose the potential well to be narrow and localized close to x = as in Fig. 1. The 
under-barrier "time" To is expressed through the particle energy E in the well 



E = 



m f dx 



+ V(0) 



(13) 



TO 
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which is weakly violated by a small non- stationary field. Certain semiclassical conditions 
(not extremely small and sufficiently slow varying S(t)) should be fulfilled. 

Eqs. (9-13) constitute the method of classical trajectories in tunneling under non- 
stationary conditions. This method relates to calculation solely of the maximum in time 
Wmax of the probability W(t) of outgoing particles. In contrast to a static barrier, the 
method of classical trajectories is not trivial in application to a nonstationary case. This 
can be seen from that the non-stationary field of the type £ cos Qt in imaginary time be- 
comes proportional to cosh Qt and may be very big. The same relates to the field (4) which 
is singular in imaginary time. This means, formally, that arbitrary weak amplitude £ of a 
nonstationary field may produce a big effect. The arguments in II also show an increase of 
an effective £ in tunneling, but it is clear physically that an effect of an extremely small 
nonstationary field is negligible and the method of classical trajectories should not work in 
this case. As one can see, the method of classical trajectories is delicate and requires an 
accurate justification. 

In the papers [25,26] for the certain particular potential well and nonstationary pulses 
the wave function has been found in the form of expansion with respect to powers of h 



where S is the classical action. The semiclassical formalism, developed in [25,26], confirms 
the method of classical trajectories and coincides with the general scheme (9-13). Neverthe- 
less, it would be extremely instructive to find any exactly solvable case of tunneling under 
nonstationary conditions in order to see how the method of classical trajectories follows not 
from a semiclassical formalism but from an exact mathematical solution. In IV this program 
is completed. 

IV. EXACT SOLUTION OF A TUNNELING PROBLEM UNDER 
NONSTATIONARY CONDITIONS 

Let us consider an one dimensional two band semiconductor in an external homogeneous 
electric field £ + £(t), where £q is the static component. The Schrodinger equation has the 
form 



where p = —ihd/dx, e g is the energy gap, and c is the velocity at big momentum. The square 
root in Eq. (15) may have two signs according to two energy bands. In a static electric field 
£o the inter-band quantum tunneling is possible which is called Zener breakdown [28]. The 
incident flux of particles from left to the right in Fig. 4 is mainly reflected back from the 
tilted energy gap but a small fraction penetrates the other band and goes to +oo with the 
probability [28] 



tp(x, t) = [a(x, t) + hb(x, t) + ...] exp — S(x, t) 



(14) 




(15) 




(16) 
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where t = a/c and 2a = e g /£ is the tunneling length in Fig. 4. Eq. (16) holds under the 
semiclassical condition 



h < e g to (17) 

The peculiarity of this problem is that it can be solved exactly with a nonstationary field by 
making the Fourier transformation with respect to x. In the new variables, momentum q and 
time t, the Schrodinger equation (15) is of the first order with respect to derivatives d/dq 
and d/dt and can be solved by the method of characteristics. If to measure the coordinate 
in units of a and time in units of t Q , the solution of Eq. (15) has the form 

^ t)= r^ exp {^it [qx - Ft{q)] } (is) 

where 



/<? / r ru-q 

duf t (u,q); f t (u,q) = Jl+ u + J dsh(s + 1) 



ru—q 

(19) 



and h(t) = £(t)/£ . Below we consider the case when the nonstationary field is much less 
compared to the static value 

£(t) < £ (20) 

If £(t) is an analytical function of t, F t (q) is an analytical function of the complex argument 
q and it has only two branch points at q — ±i. According to this, Eqs. (18-19) can be 
written through the complex paths 

^ (M)= XS exp {^ [ ^~ Ft(g)] } ; Ft{q)= Sc duft{u ' q) (2i) 

where the contour C on the complex plane of q and the contour C q on the complex plane 
of u are shown in Fig. 5. Two cuts are denoted in Fig. 5(a) by solid lines and the regular 
branch of the function F t (q) is chosen in order to get the limit 

F t (q) = - f du\f t (u,q)\; q - -oo (22) 

J — oo 

The integral (22) is divergent and has to be cut off at the lower limit by a big negative 
value which sets an irrelevant constant phase shift. Under the semiclassical condition (17) 
the g-integration in Eq. (21) goes mainly in the vicinity of the saddle point(s) determined 
by the equation 

x=^ F t (q) (23) 

Suppose the nonstationary field £{t) — > at t — > — oo when the particle energy is zero and 
the particle goes to the tilted energy gap in Fig. 4 from the left. For x to the left of the point 
—a in Fig. 4 (x < —1 in the dimensionless units used) an influence of the nonstationary 
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field under the condition (20) is small and two saddle points q = ±\/x 2 — 1 are shown in 
Fig. 6(a). They correspond to the incident and the reflected de Broglie waves. Since the 
tunneling probability is small, at x < — 1 the wave function can be shown to have the form 

#M) ~ cos (^x 2 - (24) 

In the classically allowed regions, to the left and to the right of the tilted gap, the classical 
trajectory x c i(t) (the equation, it satisfies, is written below) is weakly violated by the non- 
stationary field under the condition (20) 

xJt) = < ' 25 

excepting that x exit ^ 1 now. The energy of the incident particle is zero and the energy 
after exit from the tilted gap conserves and equals E exit = (1 — x exit )e g /2. 

For each x and t the saddle point condition (23) determines a certain q. In the absence of 
a nonstationary field the amplitude of the outgoing wave is almost stationary. Under action 
of a pulse there is an additional outgoing wave packet at t > which goes to the right and 
keeps the constant maximum amplitude (if to neglect small quantum effects of smearing) 
at the classical trajectory (25). We are interested to find this maximum amplitude of the 
outgoing packet since it defines, with the exponential accuracy, a total number of particles 
passing the barrier. With x = x c i(t) the saddle point condition (23) is satisfied by q = t 
since h(t) is small for real t. In order to reach the saddle point at positive x, one should 
deform the contour C as shown in Fig. 6(b). This could be done with no problems since 
F t (q) has only the branch point at q = % in the upper half plane. The main contribution to 
the g-integration over the bent contour C comes from the saddle point, shown in Fig. 6(b). 
The saddle point condition (23) now reads 

x exi t = 1 - 2 / dr Jo 11 U h(ir) (26) 

Jo y/l-lr + fidnhiin)] 2 

The nonstationary field is supposed to be symmetric £(t) = £(—t), as the field (4), so that 
£(ir) is real. r is a smallest positive root of the equation 
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t + drhiir) 



(27) 



and can be called the tunneling time. If to define the transition probability W as a square 
of the ratio of the maximum amplitude of the outgoing packet and the amplitude of the 
incident wave, Eq. (18) gives with the exponential accuracy 

W ~ exp {-A) (28) 

where 




d,Tih(iTi) 







(29) 
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The action A is collected due to the g-integration in the narrow vicinity of the right-hand-side 
saddle in Fig. 6(b) where q is real and positive. For such q the contour C q in Fig. 5(b) goes 
around the branch point u = ir and leads to the relation (29). Without a non- stationary- 
field Eqs. (28) and (29) give the static result (16). 

Using of the above method of saddle point means that the length 5q ~ \fh/e g to of 
the steepest descent integration should be sufficiently short compared to the typical scale 
of the problem ((r — 9) /to in dimension units) which can be estimated from the relation 
J drh{ir) ~ 1. For the pulse (4) the condition reads 

ft «m 3 # (30) 



£gto \to ) Sq 

In presence of a nonstationary field, the semiclassical condition (17) should be supplemented 
by the additional semiclassical condition (30). According to it, the nonstationary field £ 
and its duration 9 should not be too small. Note, that under the semiclassical condition 
(30) the nonstationary amplitude £ can be still smaller than the static field £q. 

In Refs. [17,18] the method of classical trajectories was developed for the present problem 
of inter-band tunneling in presence of a nonstationary field. According to this method, one 
has to minimize the classical action of a particle with the spectrum, as in (15,) and in a 
homogeneous electric field 




\dt J 

to find a classical trajectory. It satisfies the equation 

dx cl (t) t + fidsh(s) 



(^) 2 + ^ + Ht)}> (3D 



dt 

t + J dsh(s) 



(32) 



The classical trajectory starts at t — at the point x = x exit in Fig. 4. At the "moment" ir 
it reaches some point inside the gap where two branches merge and, turning on the opposite 
side of the cut, it comes at t = in the point x = —a. Within the formalism, developed in 
Refs. [17,18], in order to obtain a tunneling probability, one should substitute the classical 
trajectory into the action (31) and to integrate over the branch point ir . This procedure 
leads to the result (29). 

For the nonstationary pulse (4) under the condition (20), the tunneling time To ~ 1 at 
to < 9 and r ~ 9 /to at 9 < t . The action (29) becomes of the form 



e 



g to JV 2 ; t < 9 



A = ^< _ / 2 . . (33) 



ft \ arcsin {9 /to) + {9 /to) yj 1 - (9/t ) 2 ; 9 < t 



This expression weakly depends on £ under the semiclassical condition (30) since the role of 
the nonstationary field (4), due to its singularity in imaginary time, is only to set the new 
tunneling time r . Obtained results for inter-band tunneling under nonstationary conditions 
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may be interpreted as positive (positive energy transfer) photon-assisted tunneling shown in 
Fig. 4. In this problem cannot be such phenomenon as Euclidean resonance [26], associated 
with a negative energy transfer. 

Analogously one can consider the monochromatic £ cos Vtt [17,18] and the Gaussian 
£exp(— Q 2 t 2 ) nonstationary fields which become exponentially big in imaginary time and to 
establish how much they can grow within the formalism of classical trajectories. 

One can conclude now, that for inter-band tunneling in a nonstationary field under 
semiclassical conditions the method of classical trajectories follows from the exact analytical 
solution. The bridge between exact theory and classical trajectories enables in problems of 
tunneling to treat a nonstationary field in complex time. 



V. A NEW ENHANCED TUNNELING THROUGH STATIC BARRIERS 

Tunneling through one-dimensional static barrier is a well known problem described in 
almost every text book on quantum mechanics and the tunneling rate, with exponential 
accuracy, is given by WKB formulas (1) and (2). Tunneling through a two- or many- 
dimensional potential barrier was considered in a variety of publications [30-37]. The condi- 
tion V(r) = Eq, where E is the total energy, determines certain surfaces in the coordinate 
space. For a two-dimensional case it is shown in Fig. 7. Tunneling probability from the 
classically allowed region A to another such region B is given by Eq. (1) where the exponent 



A (A ^B) = ^ 



pro 


m 1 


/ dr 




Jo 


~2 { 



(34) 



is expressed through the classical trajectory r(r) satisfying Newton's equation in imaginary 
time 

d 2 f _ dV(r) 
dr 2 dr 

The trajectory connects different surfaces, being perpendicular to them at the connected 
points, terminates at the border of B at t = 0, and starts at the border of A at t = ir$. The 
time t is defined by the total energy E . 

The above general scheme is applicable to a calculation of tunneling rate regardless of 
dimensionality of the problem (we do not consider a possibility of caustics between A and B) . 
Nevertheless, when the dimensionality is bigger than one, another way may exist to tunnel 
through a static potential barrier. This way differs from the normal scheme and may result 
in an enhanced tunneling rate. Any tunneling process obeys quantum mechanical rules: a 
tunneling particle at each point emits partial de Broglie waves which interfere and produce 
a final outgoing wave. This wave is exponentially small for semiclassical one-dimensional 
barriers. In two- or multi-dimensional case the tunneling motion can be influenced by 
other (non-tunneling) motions leading to different conditions of interference. As a result, 
the outgoing wave may be proportional to the enhanced exponent, which is much bigger 
than normal one or even of the order of unity. The phenomenon of enhanced tunneling is 
considered below. 
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A. Tunneling in a quantum wire 



The mechanism of enhanced tunneling through a static barrier can be shown in the 
case of a long and narrow quantum wire when the radial motion (x-direction) is strongly 
quantized but the motion along the axis (y-direction) is almost classical. We do nod discuss 
a detailed link to real quantum wires since the goal is solely to demonstrate a new mechanism 
of tunneling through static barriers. To model this situation, suppose the particle to move in 
the two-dimensional potential V x (x) + V y (y) plotted in Fig. 8. The particle can tunnel from 
the energy level E in the radial direction. In two dimensions the potential energy is shown 
in Fig. 9. The particle tunnel from the region of small x to the region of big x. When the 
total potential energy is simply a sum of ^-dependent and y-dependent potentials, variables 
x and y are separated and tunneling along re-direction occurs as in a one-dimensional case. 

B. Mapping of the stationary problem onto dynamical pulses 

If to add some "interaction" potential V int (x,y) to the system, the total potential 

V(x,y) = V x (x)+V y (y)+V int (x,y) (36) 

is not reduced to a sum of x- and y-term and motions in the x- and y-direction are no 
more independent. The both barriers in Fig. 8 are weakly penetrated and, in the first 
approximation, the particle occupies the region of x ~ Xq and positive y in Fig. 9 with 
incident and reflected fluxes along y-direction as shown in Fig. 8(b). This is equivalent 
to the problem of two one-dimensional particles with the interaction V int (x,y) when one 
particle in a well is acted by a steady flux of incident particles. This is analogous to the 
situation in nuclear physics when an incident steady flux of protons acts on the alpha particle 
occupied some energy level inside the nuclear potential well [38,39]. This nuclear problem 
of excitation of alpha particle can be considered by two methods: (i) by a solution of the 
Schrodinger equation for alpha particle and proton in the static potential or (ii) by a solution 
of the Schrodinger equation for alpha particle in a nonstationary potential created by the 
classically moved proton which reflects by nuclear Coulomb forces. The both ways lead to the 
same small exponent in the semiclassical probability. The second method is more convenient 
since the excitation of a higher level of alpha particle can be considered as one driven by 
a nonstationary force, acting during a finite time [9]. So, with exponential accuracy, the 
steady problem of decay of a metastable state under an incident flux can be mapped onto 
dynamical pulses. 

Under conditions, when the motion in the y-direction is weakly violated by the motion 
in the x-direction, the total process can be considered as one-dimensional tunneling in the 
nonstationary potential V x (x) + V int [x,y(t)], where y(t) is a classical trajectory in the static 
potential V y (y), shown in Fig. 8(b). According to this, besides normal tunneling "0" in 
Fig. 8(a), corresponding to absence of the pulse, there is also the process "1" related (in 
the language of dynamical pulses) to the under-barrier photon-assisted tunneling which is 
similar to Fig. 4. Since a photon-assistance increases tunneling probability, the process "1" 
relates to the phenomenon of enhanced tunneling. In the language of static potential, one 
can say that tunneling along the radius is assisted by the motion along the axis. 
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This mapping is mentioned only for interpretation of enhanced tunneling "1" by means 
the language of dynamical pulses. In VC enhanced tunneling is considered on the basis of 
a static barrier. 



C. Probability of the enhanced tunneling 



The condition V(x,y) = E gives the curve shown in Fig. 10. To calculate a probability 
of enhanced tunneling "1", one should account the fact that in semiclassical limit the phase 
S/h of the wave function ip(x,y) ~ exp[iS(x,y)/h] is big. The function S(x,y) can be 
found from Scrodinger equation as the main term in the expansion of the type (14). For 
normal tunneling "0" , S(x, y) coincides with the classical action, varies smoothly between 
the regions A and B, and its extreme value relates to the classical trajectory denoted by the 
dashed curve in Fig. 10. For enhanced tunneling "1" the behavior of S(x,y) is much more 
complicated, since there are reconnection among different branches of S. On each branch the 
function S(x,y) coincides with some classical action according to findings of Refs. [25,26]. 

Fortunately, to find a tunneling probability it is not necessary to know ip(x,y) at all x 
and y. It is sufficient to connect phases of the wave function between certain points at the 
borders of A and B using some formal method. Such formal method for the process "1" 
(enhanced tunneling) consists of two steps: a connection A — > / and a connection / — > B 
(Fig. 10). For the first step, one can consider the vertical line at x = xq connecting the 
border of A and some point / in Fig. 10, where ~ exp(ia/h) is expressed through the 
certain branch a of the classical action, which satisfies the equation of Hamilton- Jacobi 



At x > Xq the potentials V x (x), V y (y), and Vi nt (x,y) are smooth. Since the potential V x (x) 
is abrupt at x = x , existence of the quantum level E in the well results in the following 
semiclassical condition inside the barrier in Fig. 8(a), when x tends to x from the right, 



The total energy is E = E + s, where e is the energy of the classical motion along the well. 
The y-dependence of the action at X — Xq Cclll be easily found from Eqs. (37) and (38) in the 
WKB form of the one-dimensional motion along the well. The phase at the point / (relative 
to A) in Fig. 10 is 



where yo an d yi are determined by zero of the square root in Eq. (39) at / and A respectively. 
The position of the point /, defined by the conditions x = xq and y = yo, is not flexible 
and depends only on the energy e, since the quantum level E is strictly determined by the 
potential in Fig. 8(a). The same relates to the phase a /hat the point / (relative to A). The 
value of a, given by Eq. (39), is extreme with respect to variation of positions of the ends 
of the line x = Xq in Fig. 10. 




(37) 




(38) 




(39) 
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As the second step, the point / can be connected to the region B by the classical 
trajectory which provides an extreme of the phase of B (relative to /) and is based on the 
expression A (f — > B), as Eq. (34), but the trajectory connects now the regions / and B. 
The rate of enhanced tunneling "1" exp(— Ai) is given by the total phase difference and 
reads 

2? a 

A, = A (f -+B)-— (40) 

where a is defined by Eq. (39). As soon as normal tunneling is described solely by a classical 
trajectory (34)- (35), enhanced tunneling, according to Eq. (40), requires more complicated 
description. The total tunneling probability should be estimated as 

W ~ max{exp(-A ); exp(-Ai)} (41) 

To calculate A (normal) and Ai (enhanced) one should specify the potential energy in the 
form (36) and use Eqs. (34), (35), (39), and (40). 

The second term in Eq. (40), which is real and negative, works versus the big first term 
representing the double tunneling through the barriers in Figs. 8(a) and 8(b). The second 
term in Eq. (40) can play a crucial role in reduction of A\. Enhanced tunneling "1" may 
be interpreted as a cooperative motion along a trajectory in the x-direction and along a 
bounce (forth and back) in the y-direction. One should emphasize again, that the above 
semiclassical method gives only the formal calculation of the tunneling probability with 
exponential accuracy, but it does not relate to a real behavior of the wave function if)(x,y). 
For example, the formal density at the point / can be big, the real wave function drops 
exponentially in that region. 



VI. DISCUSSION AND CONCLUSIONS 

There are two issues in this paper. 

The first one is that the problem of inter-band tunneling in a semiconductor (Zener 
breakdown) in a nonstationary and homogeneous electric field is solved exactly. On the basis 
of the exact analytical solution the formation of an approximation of classical trajectories 
is studied. This semiclassical formalism is very effective in tunneling under nonstationary 
conditions since it allows to reduce the problem to a solution of the Newton equation in 
complex time. The semiclassical approach, obtained from the exact solution, is of the same 
type as in other problems of nonstationary tunneling which do not allow exact solutions, 
for example, decay of a metastable state, penetration through a potential barrier, and over- 
barrier reflection. 

The second issue relates to different types of tunneling through a static potential barrier 
if dimensionality of the problem is bigger than one. A mult i- dimensional case has an essential 
feature which differs it from one-dimensional problem. Namely, a tunneling coordinate (x) 
can be influenced by non-tunneling one (y) resulting in a different condition of interference of 
emitted partial waves during tunneling through the potential V(x,y). A mult i- dimensional 
tunneling can be interpreted as one-dimensional but in the nonstationary potential V[x, y(t)], 
where y(t) is the classical non-tunneling motion. This mapping onto dynamical problem 
enables to establish different types of tunneling. 
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As well known, there is normal tunneling through a mult i- dimensional barrier related sim- 
ply to a classical trajectory in imaginary time and connected two classically allowed regions. 
This mechanism was always employed by authors in consideration of mult i- dimensional tun- 
neling. Normal tunneling corresponds to the case when the tunneling particle does not 
absorb quanta of the nonstationary field. 

There is also a process when particle absorbs quanta (photons) from the nonstationary 
potential V[x, y(t)]. This is a photon-assisted tunneling with an enhanced probability. In 
the language of the static potential barrier, this process corresponds to enhanced tunneling. 
This is a new mechanism of quantum tunneling through a static multi-dimensional bar- 
rier. Depending on shape, sign, and value of the potential V(x,y), two types of enhanced 
tunneling are possible: (i) with absorption of energy by the tunneling motion (positive as- 
sistance) or (ii) with emisssion of energy (negative assistance). The positive assistance was 
considered in Ref. [25] and in this paper for Zener breakdown. The negative assistance 
was considered in Ref. [26]. It is remarkable that the method of calculation of probability 
of enhanced tunneling, proposed in this paper, is applicable equally to both positive and 
negative assistance. 

As it was found in Ref. [26] for tunneling through a nonstationary barrier, the enhanced 
tunneling probability may be not exponentially small under conditions of negative assistance. 
This phenomenon is called Euclidean resonance which results from influence of nonstation- 
ary field on interference of emitted partial waves during tunneling. Euclidean resonance 
can also take place in tunneling through a static non-one-dimensional barrier. To get this 
phenomenon, one should adapt barrier parameters and the particle energy to reach the con- 
dition A\ = in Eq. (40), as in the case of a nonstationary barrier [26]. The above theory 
is valid if exp(— Ai) <C 1 and when this condition breaks down one should use a formalism 
generic with the multi-instanton approach. Nevertheless, the presented theory, as any other 
approach, can be approximately extended up to the limit of its validity exp(— A\) ~ 1 wich 
relates to Euclidean resonance. 

Enhanced tunneling may result in a dramatic increase of a tunneling rate. For example, 
for an almost classical barrier, according to the normal mechanism, the tunneling probability 
is calculated to be 10~ 100 , but enhanced tunneling, under conditions of Euclidean resonance, 
results in the probability of 10 -2 . Euclidean resonance is a phase phenomenon and any 
effect of dephasing (friction, for example) may destroy it. To avoid this destruction, the 
attenuation time due to friction should be bigger compared to the under-barrier time. 

In conclusion, a new mechanism of enhanced tunneling through static non-one- 
dimensional barriers is established in addition to well known normal tunneling. Under 
certain conditions on shape of an almost classical barrier and particle energy, the probabil- 
ity of quantum tunneling through it may be not exponentially small. 
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FIGURES 




FIG. 1. The path of tunneling is denoted by the dashed line, xt is the classical turning point. 
E is the energy of the metastable state, V is the barrier height, and a is the typical potential 
length. 
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V(x) 




FIG. 2. (a) The particle can absorb a quantum and tunnel in a more transparent part of the 
barrier with the energy E + hft; (b) the process of the multi-quanta absorption with the subsequent 
tunneling. 
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wave function 




FIG. 3. The outgoing wave function at the turning point xy has a maximum amplitude under 
the influence of the external pulse. 



energy 




FIG. 4. An inter-band tunneling (Zener breakdown) in a semiconductor with the energy gap 
e g in a static electric field So, shown by the open arrow. The tunneling length is 2a = e g /£o. With 
a nonstationary field the path inside the forbidden region is shown by the full arrow. 
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FIG. 5. (a) The cuts are shown by solid lines; (b) the contour C q terminates at the point q. 
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saddle 
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FIG. 6. (a) At x < —a two saddle point are involved which correspond to the incident and the 
reflected de Broglie waves; (b) for positive x the saddle point relates to a particle after tunneling. 
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X 



FIG. 7. Thick curves relate to the condition V(x,y) = Eq. Normal tunneling between two 
classically allowed region A and B, where V(x,y) < Eq, is described by the classical trajectory in 
imaginary time (thin curve) which meets the borders of A and B normally. 




FIG. 8. The potential energy of a particle in a quantum wire V x (x) + V y (y) corresponds to the 
quantized motion in the radial x-direction (a) and an almost classical motion in the axis y-direction 
(b). 
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y 

FIG. 9. The potential energy V(x,y) of a particle in the quantum wire in two dimensions. 



y 




FIG. 10. Thick curves are determined by the condition V(x,y) = Eq. Normal tunneling "0" 
simply relates to the classical trajectory denoted by the dashed curve. Enhanced tunneling "1" 
is more complicated and consists of two steps: A —>■ f and the classical trajectory / —>■ B. The 
vertical line goes from y = yi down to y = yo- 
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